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PREFACE 


This  report  was  prepared  by  the  Air  Force  Engineering  and  Services 
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Faculty  Research  Program. 

The  report  covers  work  performed  between  21  May  and  29  July  1984.  The 
AFESC/RDCS  Project  Officer  was  Capt  Paul  L.  Rosengren,  Jr. 

This  study  shows  that  Finite  Element  analysis  can  be  a  useful  tool  in 
examining  the  validity  of  the  results  of  centrifugal  model  testing  as  long  as 
the  constitutive  relationships  for  soil  are  representative  of  actual  behavior. 


This  report  has  been  reviewed  by  the  Public  Affairs  Office  (PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NT1S  it 
will  be  available  to  the  general  public, including  foreign  nationals. 
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SECTION  I 
INTRODUCTION 


A.  PURPOSE 

Soil-structure  interaction  problems  have  lone  been  of  interest  and 
concern  for  geotechnical  engineers.  Although  rapid  progress  has  been 
made  in  solving  these  problems,  their  application  to  concrete  culverts 
embedded  in  embankments  has  been  slow.  Due  to  an  increased  number  of 
concrete  culvert  projects  in  industry,  state  highway  departments,  and 
federal  agencies  in  the  United  States  during  recent  years,  research 
could  improve  the  design/analysis  of  concrete  culverts  and  result  in 
better  structural  stability,  and  time  and  cost  savings. 

B.  OBJECTIVES 

One  currently  used  design  method  for  buried  culverts,  originally 
developed  by  Marston  (Reference  1)  early  in  the  century,  has  been 
modified  by  Spangler  (Reference  2)  and  Costes  (Reference  3).  According 
to  this  method,  the  design  of  buried  culverts  in  a  soil  fill  depends  on 
two  empirical  factors,  settlement  ratio  and  load  factor.  However, 
primary  shortcomings  of  this  method  are:  (1)  the  settlement  ratio  is  not 
known  before  the  culvert  is  built;  and  (2)  the  designer  must  predict 
this  parameter  from  a  few  past  experimental  observations  or  by 
engineering  judgment.  Therefore,  there  is  a  need  for  more  rational 
design  methods  on  a  sounder  theoretical  basis. 

Perhaps  the  most  ideal  approach  for  obtaining  information  for  the 
development  of  new  design  methods  is  full-scale  model  testing.  A  full- 
scale  model  with  the  necessary  instrumentation  (i.e.,  soil  stress 
meters,  settlement  gauges,  and  strain  gauges,  etc.)  could  give  the  best 
results  for  estimating  prototype  behavior  of  buried  culverts. 

Unfortunately,  full-scale  model  testing  has  serious  najor 
drawbacks:  mainly,  cost  and  time  of  construction  and  operation.  Because 
of  these  reasons,  the  centrifugal  modeling  technique  is  becoming  a 
favorite  testing  method  in  geotechnical  engineering.  This  technique 
offers  a  comparatively  inexpensive  and  easy  way  of  obtaining  essentially 


the  saine  data  as  could  be  obtained  from  field  tests.  It  also  provides 
opportunities  for  studying  the  effect  of  individual  parameters,  and 
collecting  more  data  than  field  tests  usually  yield. 

However,  at  present,  there  is  limited  capability  to  add  or  remove 
soil  in  centrifugal  flight.  Thus  a  typical  test  of  a  model  embankment 
involves  subjecting  that  embankment  to  a  gradually  increasing 
centrifugal  field,  while  maintaining  it  at  its  final  shape.  Although  the 


models  can  be  made  large  enough  to  include  blocks  of  undisturbed  and 
compacted  soil  which  are  representative  of  field  conditions  (References 
A),  an  objection  could  be  raised  that  the  stress  paths  traversed  by  soil 
elements  during  the  test  are  unrepresentative  of  those  \diich  apply  in 
the  field. 

Therefore,  this  necessitates  computational  methods  which  examine 
the  validity  of  such  objections.  Finite  element  techniques  can  be  used 
to  create  mathematical  models  which  model  both  the  laboratory  and  the 
field  situations,  and  allow  the  two  to  be  compared.  As  long  as 
constitutive  equations  for  the  soil  used  in  the  mathematical  model  are 
representative  of  real  ones,  valid  comparisons  should  be  possible. 

Therefore,  the  objectives  of  the  study  ares 

1.  to  investigate  the  influence  of  the  constitutive 
equations  for  soil  on  the  culvert  response  in 
finite  element  analysis,  and 

2.  to  compare  the  results  of  the  finite  element  analysis 
to  the  results  of  the  centrifugal  model  testing  of  the 
concrete  culvert. 


SECTION  II 

CENTRIFUGAL  MODEL  STUDY 


The  centrifugal  model  testing  of  the  buried  concrete  culvert  in 
sand  was  performed  by  James  and  Larsen  (Reference  5)  at  Cambridge 
University,  and  only  brief  descriptions  for  the  model  testing  procedures 
are  summarized  here. 

A  strong  17.8- by  15.2- by  5.9-inch  rectangular  box  was  mounted  on 
the  centrifuge.  The  box  consisted  of  a  frame  made  of  hollow  sections 
welded  together,  the  frame  was  bolted  to  a  back  plate  of  mild  steel  and 
the  front  of  the  box  was  made  of  a  Plexigla^ plate,  which  permitted 
observation  of  the  model. 

After  mounting  the  strong  box,  the  pipe,  which  was  cast  of 
microconcrete  (i.e.,  concrete  with  scaled-down  aggregates),  was  placed 
1.2  inches  away  from  the  bottom  of  the  box.  The  load  cells  were  mounted 
in  the  pipe  wall  to  measure  the  normal  tractions,  and  strain  gauges  were 
glued  on  the  inside  and  outside  of  the  pipe  wall  to  measure  the  bending 
moments  and  thrusts  of  the  pipe.  Sand  was  poured  from  a  hopper  into  the 
box  with  approximately  constant  height  of  fall  and  rate  of  pouring.  To 
reduce  side  friction  in  the  box, glass  plates  were  used  between  the  sand 
and  steel  plate  and  between  the  sand  and  the  Plexiglas?  plate. 

The  speed  of  the  centrifuge  was  then  raised  in  steps  of  5  g 
increments  and  the  final  load  cell  readings  and  strain  gauge  readings 
were  taken  at  35  g. 

Figure  1  shows  the  geometry  and  boundary  conditions  of  the  model. 


3 


Figure  1.  Geometry  and  Boundary  Conditions 
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SECTION  III 

FINITE  ELEMENT  PROGRAM 


The  finite  element  program  used  in  this  study,  CANDE  (Culvert 
ANalvsis  and  DEsign)  was  developed  by  Katona  (References  6  and  7). 
Although  CANDE  can  be  operated  in  a  design  or  analysis  mode,  only  the 
analysis  mode  was  utilized  in  this  study.  The  main  features  of  CANDE, 
which  are  used  in  the  analysis,  will  be  briefly  described  in  this 
report . 

The  assumptions  common  to  all  solution  levels  in  CANDE  are:  (1) 
Two-dimensional  analysis  with  plane  strain;  and  (2)  small  strain  theory. 
There  are  three  solution  levels  in  CANDE.  Level  1  corresponds  to  the 
elastic,  closed  solution  while  Levels  2  and  3  are  both  finite  element 
(FE)  solutions.  The  difference  between  Levels  2  and  3  is  in  the  way  the 
FE  mesh  is  defined.  Level  2  generates  automatic  and  symmetric  meshes  for 
a  circular,  elliptical,  or  box  culverts,  whereas  Level  3  allows  the  user 
to  define  a  mesh.  Since  the  model  culvert  used  in  this  study  was 
symmetric  and  circular,  Level  2  was  used. 

Five  material  characterizations  are  available  for  various  pipe 
materials  in  CANDE.  These  materials  are:  (1)  steel;  (2)  aluminum;  (3) 
plastic;  and  (4)  concrete.  In  this  study  the  concrete  material 
characterization  with  a  bilinear  stress-strain  relationship  was  used  to 
represent  the  concrete  material  properties. 

Three  basic  elements  are  contained  in  the  program:  a  quadrilateral 
element  for  representation  of  the  soil,  an  interface  element  for  slip  of 
elements,  and  a  beam-column  element  for  the  culvert. 

For  representation  of  soil  behavior,  four  types  of  constitutive 
relationships  are  incorporated  witli  the  program.  These  are:  (1)  linear 
elastic;  (2)  overburden-dependent;  (3)  extended  Hardin’s  model;  and  (4) 
Duncan’s  hyperbolic  model.  To  investigate  the  influence  of  the 
constitutive  relationship  for  soil  on  the  culvert  response  in  the 
numerical  analysis,  both  a  linear  elastic  model  and  a  nonlinear  model 
(Duncan’s  hyperbolic  model)  for  soil  were  used  in  the  analyses.  The 
results  for  each  representation  are  presented  in  this  report.  In 


Reference  5  only  dense  sand  with  a  void  ratio  of  0.52  was  specified, 
soil  properties  for  the  elastic  (FE)  soil  model  were  assumed  as:  Young's 
modulus  of  7,000  psi  and  Poisson's  ratio  of  0.3.  A  "SW"  soil  type  at  a 
compaction  of  level  of  95  percent  in  Table  1  was  used  for  Duncan's  soil 
model.  These  values  are  within  an  expected  range  as  recommended  by  Lambe 
(Reference  8)  and  Duncan  (Reference  9),  respectively. 

The  finite  element  grid  with  boundary  conditions  used  in  the 
comparison  analysis  is  shown  in  Figure  2.  The  culvert  was  represented  by 
10  beam-column  elements,  the  soil  was  represented  by  86  quadrilateral 
elements,  and  the  slip  model  was  represented  by  11  interface  elements. 


TABLE  1.  REPRESENTATIVE  PARAMETER  VALUES  OF  THE  MODIFIED  DUNCAN  MODEL  (REFERENCE  9) 
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RC  =  relative  compaction,  in  percent 


Number  of  Nodes  =132 
Number  of  Beam  Elements  =  10 
Number  of  Quadrilateral  Elements  =  86 
Number  of  Interface  Elements  =  11 
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Figure  2.  Symmetric  Mesh 


SECTION  IV 
RESULTS 


The  quasi-theoretical  results  obtained  from  the  finite  element 
code,  CANDE,  are  compared  to  the  results  of  the  centrifugal  testing 
model.  The  compared  quantities  are i  (1)  normal  traction  on  the  culvert 
periphery;  (2)  peripheral  distribution  of  bending  moment;  and  (3) 
peripheral  distribution  of  thrust. 

Figure  3  shows  graphical  comparisons  between  finite  element  results 
and  centrifugal  model  testing  observations  for  normal  tractions  acting 
at  the  culvert  periphery.  Excellent  comparisons  are  observed  from  the 
results  of  numerical  predictions  with  Duncan's  soil  model  and  the 
results  of  the  centrifugal  model  testing.  Trends  of  shape  changes  are 
similar  and  the  magnitudes  are  in  excellent  agreement.  In  particular, 

0  0  O 

the  magnitudes  around  the  positions  of  30  ,  90  (springing  line),  144  , 
and  162  from  the  crown  are  almost  identical  to  each  other. 

Yet,  discrepancies  exist  between  the  numerical  predictions  with  the 
linear  soil  model  and  other  results  obtained  from  the  centrifugal 
testing  model  and  numerical  predictions  with  nonlinear  soil  model  -  even 
though  the  trend  of  shape  change  is  similar. 

Figure  4  shows  internal,  peripheral  bending  moment  distributions  of 
the  culvert.  As  shown  in  from  the  figure,  moderately  good  comparisons 
are  observed  between  the  results  of  numerical  analyses  and  the  results 
of  the  centrifugal  testing  model.  Shape  changes  are  similar  and  the 
magnitudes  are  in  good  agreement  except  around  the  invert.  Positive 
moment  in  the  figure  is  that  which  produces  tension  at  the  inner  fibers 
of  the  concrete. 

Figure  5  compares  the  results  of  finite  element  analyses  and  the 
centrifugal  model  study  of  circumferential  thrusts  of  the  culvert. 
Numerical  predictions  and  model  measurements  are  similar  in  shape,  but 
do  not  agree  well  in  magnitudes  around  the  crown.  Correlations  are  very 
good  around  the  springing  line  and  the  invert. 

Again,  the  results  of  CANDE  with  an  elastic  soil  model  are  somewhat 
conservative  when  compared  to  other  results. 
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Figure  5.  Peripheral  Distribution  of  Thrust 


SECTION  V 
RECOMMENDATIONS 


This  study  evaluates  the  capabilities  of  a  finite  element  method 
for  analyzing  the  performance  of  buried  structures  by  comparing  the 
predicted  and  measured  behavior  of  a  concrete  culvert.  The  findings  and 
recommendations  for  further  research  based  on  this  study  arc  as 
follows: 

1.  The  results  of  the  finite  element  analyses  compared  well 

in  shape  with  the  results  of  the  centrifugal  testing  model. 
Analysis  with  Duncan's  hyperbolic  relationships  for  soil 
exhibits  a  more  realistic  approach  for  predicting  soil  and 
culvert  responses;  the  results  were  much  closer  to  the  results 
of  the  centrifugal  testing  model  in  shapes  and  magnitudes. 

2.  Significant  economic  savings  (time  and  costs  of  construction 
and  operations)  could  be  obtained  by  using  the  centrifugal 
modeling  technique  to  study  the  behavior  of  prototype  buried 
culverts  if  results  of  the  centrifuge  model  testing  agree 
well  with  field  measurements. 

3.  Additional  economic  savings  could  be  obtained  when  numerical 
(finite  element)  analyses  are  used  with  realistic  constitutive 
models  for  soils  and  culverts. 

A.  Additional  studies  (i.e.,  box,  arch,  or  elliptical  culverts 
with  different  bedding  and  boundary  conditions)  should  be 
conducted  for  further  verifications  before  general  acceptance 
can  be  awarded. 
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